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Abstract 

o ; 

It is demonstrated using Monte Carlo simulation that in different nucleus— nucleus 
Ph. collision samples, the increase of the fluctuation of event factorial moments with de- 

creasing phase space scale, called erraticity, is still dominated by the statistical fluctu- 
ations. This result does not depend on the Monte Carlo models. Nor does it depend on 
^ ■ the concrete conditions, e.g. the collision energy, the mass of colliding nuclei, the cut of 

phase space, etc.. This means that the erraticity method is sensitive to the appearance 
of novel physics in the central collisions of heavy nuclei. 
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It is generally believed that through the collision of heavy nuclei at ultra-high ener- 
gies big systems with very high energy density fl[] might be produced. In these systems 
novel phenomena, such as colour deconfmement 0, chiral-symmetry restoration ||, 
discrete-symmetry spontaneous-breaking [|J, etc., are expected to be present and dif- 
ferent events might be governed by different dynamics. With this goal in mind, the 
event-by-event (E-by-E) study of high energy collisions has attracted more and more 
attention ||. 

A well known example of E-by-E fluctuation is the dynamics of self-similar cascade, 
which results in a fractal system, and the dynamical probability-distribtion fluctuates 
E-by-E U . Such kind of self-similar dynamical fluctuations can be studied by means 
of the method of normalized factorial moments (nfm) |J. The latter are defined as 

jp rn r\ 1 V" ( n m(n m - 1) • ■ ■ (n m - q + 1)) 

= m\ ' {) 

m=l x 

where a region A in 1-, 2- or 3-dimensional phase space is divided into M cells, n m is the 
multiplicity in the mth cell, and (• • •} denotes vertically averaging over the event sample, 

(•••) = ^E(-). (2) 

i=l 

AT is the number of events in the sample. If self-similar dynamical fluctuations exist, the 
nfm will possess an anomalous scaling property with the diminishing of phase space scale (or 
increasing of partition number M), 

Fg(M) cx (M)*« (M — > oo) . (3) 

Recently the predicted anomalous scaling of nfm, Eq.(3), has been successfully observed in 
experiments (For a review see ||). 

In Eq.(l) the vertical average (• • •) over the event sample precedes the horizontal average 
(1/M) J2m=i(' ' ') over the M bins. The NFM defined in this way is sometimes refered to as 
vertically averaged factorial moment and denoted by Fq V \M). 

f(JO = Jg t < "" ( ""-' ) . " +1>> . (4) 
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Alternatively, one can also reverse the order of the two average processes, i.e. doing the 
horizontal average first, and define horizontally averaged factorial moment as 



It can be shown that if the vertical nfm has the anomalous scaling property, Eq.(3), then the 
horizontal NFM will have the same property. 

Note that in the definition Eq.(5) of horizontal NFM an average over the event sample has 
been made for the event normalized factorial moment Fq C ^ (M) (efm) defined as 

F^{M)- m , (6) 



where n m is the multiplicity in the mth cell of that event. Therefore, it is natural to ask the 
question: How about the E-by-E fluctuation of efm Fq^l 

Cao and Hwa |10| propose to quantify this fluctuation by the normalized moments 



C P , q = ((^y), ^=F^/(F^) (7) 

of F 9 (e) . If C p>g has a power law behavior as the division number M goes to infinity 

C M (M) oc M^ q M -> oo, (8) 

then the phenomenon is referred to as erraticity, and is characterized by the slope fi q of ip q (p) 
at p = 1 

d_ 

dp 



(9) 

p=l 



which is called entropy index. Define 

q dp 



= {^ln^), (10) 

P =i 



then the entropy index [i q can be calculated through 

The usefulness of erraticity, or entropy index, in the study of E-by-E fluctuation is limited 
by the fact that this behaviour is dominated by statistical fluctuations when the multiplicity 



is low [11]. Only for high multiplicity events, as for example in the central collisions of heavy 
nuclei, the "entropy index" coming from statistical fluctuations becomes very small and the 
dynamical effect can be expected to show up 

In the present letter this problem is studied in some detail using the Monte Carlo gen- 
erators Fritiof and Venus. It will be shown that within the framework of these models the 
statistical fluctuations still dominate the erraticity behaviour of central nuclear collisions, 
even though the multiplicity is as high as several hundreds to several thousands. What is 
interesting is that this dominance of statistical fluctuations does not depend on the model 
used. Neither does it depend on any physical condition, e.g. the collision energy, the mass 
of the colliding nuclei, the cut of phase space, etc. This means that the erraticity method 
has the peculiar property that it is able to filter out all the concrete physical conditions used 



in data analysis and therefore may be used as a sensitive signal for the appearance of novel 
physics. 

We start from the study of Pb-Pb collisions. Two samples are generated using Fritiof for 
the incident energies 158 and 500 A GeV, each consisting of 10 000 events. The phase space 
regions used for the study of erraticity behaviour are listed in the first 3 rows of Table I. The 
collisions are central in the sense that the impact parameters lie between and 0.5 fm. 



Table I The phase space region, average multiplicity (N) and entropy 
index /i2 in Fritiof Monte Carlo of Pb-Pb collisions 





Incident energy (A GeV) 


158 


500 


y 


[1,2] 


[0,1] 


[0,2] 


[-2,2] 


[-2,2] 


Pt (GeV/c) 


[0,10] 


[0,10] 


[0,10] 


[0,10] 


[0,10] 


<P 


[-7T,7r] 


[-7T,7r] 


[— 7T, 7r] 


[-7r,vr] 


[-7T,7r] 


(N) 


286.1 


407.2 


693.2 


1397.9 


1677.7 


f*2 


0.487 


0.273 


0.0857 


0.0167 


0.00856 



In order to elliminate the effect of non-flat average distribution, the phase space variables 
y, pt, (p are transformed into the corresponding cumulant forms |l3| X y , X pt: X v as usual. 
After the transformation, the phase space regions of all three X a (a = y,pt, f) become [0,1]. 

In calculating the efm, the phase space region in each direction is divided into M sub- 
cells. The total number of sub-cells in the 3-D phase space region A is M3D = M 3 . The 
log-log plots of the event-space moment C Pi 2 of efm versus M3D are shown in the left column 
of Fig's.l and 2 for p = 0.5, 0.7, 0.9, 1.0, 1.1, 1.5, 2.0, respectively. 

The derivatives £2 of C Pt 2 at p = 1 versus log M3D are plotted in the right column of 
Fig's.l and 2. The entropy indices [i<i are then obtained as the slope of £2 versus log M3D at 
large M. The results are listed in the last row of Table I. 

It can be seen from the figures that the log-log plots of C p .2 versus M3D have similar 
shape for all the cases but only with different scales. This means that erraticity exists in all 
the cases with different strength, characterized by the different values of entropy index /i. A 
regularity that can easily be observed from Table I is that the entropy index [i decreases with 
increasing average mutiplicity (N). 

The dependence of /i2 on (N) is plotted in Fig. 3. The full line in this figure is the result 
of pure statistical fluctuations taken from Ref. [12]. Our results lie well above this line, which 
seems to indicate that some dynamical effect shows up. However, this conclusion cannot 
be drawn because the full line was obtained from the pure-statistical-fluctuation model in 
one-dimensional phase space [p.2|], while our results are for 3-dimensional case. 
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Fig.l Log C P) 2 and £2 versus log M for Pb-Pb collisions at 158 A GeV obtained by the Fritiof 
generator. The rapidity regions (in cm. s.) in (a),(b),(c),(d) are: y £ [—2,2], [0,2], [0,1], [1,2], 
respectively. The transverse momentum region is p t € [0, lOGeV/c] and the azimuthal region 

is [if € — 7T, 7r] 



In order to make a faithful comparison between the results from the Fritiof generator 
and the pure-statistical-fluctuation case, we construct models of pure statistical fluctuations 
in 1-, 2- and 3-dimensions, respectively. For illustration, consider the 2-D model. Let X a 
and Xf, denote the two (cumulant) variables. For each particle in an event take two random 
numbers distributed uniformly in the region [0,1] as the values of X a and Xj, of this particle. 
Repeating N times, the X a and Xj, values of all the N particles in the event are determined 
and a Monte Carlo event, containing only statistical fluctuations, is obtained. Constructing 
in this way J\f events, the C p „ and £„ can be calculated. Note that, by construction, for the 
characterization of each particle in the 1-, 2-, 3-D models we need 1, 2, 3 random numbers, 
respectively. Therefore, the "degree of randomness" is higher and the entropy index fi q should 
be larger for the 3-D (2-D) model than for the 2-D (1-D) ones. 

The results of the calculation shown in Fig.3 as full (1-D), dashed (2-D) and dotted (3-D) 
lines confirm the expectation. A striking fact which can be seen from the figure is that the 
results of the Fritiof Monte Carlo for Pb-Pb collisions at 158 and 500 A GeV all lie on the 
dotted line, which means that the erraticity phenomena observed in the Fritiof-Monte-Carlo 
simulation of Pb-Pb collisions at these two energies are dominated by statistical fluctuations, 
inspite of the high multiplicities. 




Fig.2 The same as Fig.l, but at incident energy 500 A GeV. 
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Fig.3 The dependence of log \ii on (N). Full circles are from Fritiof Monte Carlo. 

Full stars are from Gaussian-alpha model. Full, dashed and dotted lines are 
the results of pure statistical fluctuations in 1-, 2- and 3-D, respectively. 

In order to check whether this conclusion depends on the projectile and target nuclei 
and/or on the event generator used, similar analysis is carried out for various colliding systems 
at different incident energies using both Fritiof and Venus event generators. 

Table II The average multiplicity and entropy index of nuclear 

collisions obtained from Fritiof-Monte-Carlo for different 
projectile-targets, incident energies, rapidity regions and particle types 



colliding 


± -'mc 


rapidity 


particle 


average 


entropy 


nuclei 


(A GeV) 


region 


type 


multiplicity 


index \ii 


O-Au 


200 


[-1,1] 


charged 


104.1 


0.908 


S-Au 


200 


[-1,1] 


charged 


152.4 


0.825 


S-S 


158 


[0,2] 


charged 


96.3 


0.908 


S-S 


158 


[-2,2] 


charged 


192.5 


0.718 


Pb-Pb 


158 


[1,2] 


charged 


286.1 


0.336 


Ag-Ag 


158 


[0,2] 


charged 


360.2 


0.365 


Pb-Pb 


158 


[0,1] 


charged 


407.1 


0.236 


Pb-Pb 


158 


[0,2] 


charged 


693.2 


0.0876 


Ag-Ag 


158 


[-2,2] 


charged 


721.4 


0.0891 


Pb-Pb 


500 


[0,3] 


charged 


1069.9 


0.0338 


Pb-Pb 


158 


[-2,2] 


charged 


1397.9 


0.0196 


Pb-Pb 


500 


[-3,3] 


charged 


2169.2 


0.0071 
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Fig. 4 The dependence of [i<i on (N) from Fritiof and Venus Monte Carlo 
compared with the 3-D pure-statistical-fluctuation model. 
The phase space regions used are listed in Tables II and III. 
Full stars are from Gaussian-alpha model. 

Table III The average multiplicity and entropy index of nuclear 

collisions obtained from Venus-Monte-Carlo for different 
projectile-targets, incident energies, rapidity regions and particle types 



colliding 


^inc 


rapidity 


particle 


average 


entropy 


nuclei 


(A GeV) 


region 


type 


multiplicity 


index \ii 


H-H 


650 


[-4,4] 


all 


14 


1.8499 


Pb-Pb 


158 


[0,1] 


negative 


21 


1.509 


Pb-Pb 


158 


[0,2] 


negative 


23 


1.507 


Pb-Pb 


200 


[1,2] 


all 


26 


1.519 


O-Au 


200 


[-1,1] 


negative 


57 


1.277 


S-Au 


200 


[-1,1] 


negative 


80 


1.122 


Pb-Pb 


200 


[0,1] 


all 


154 


0.8787 


Pb-Pb 


200 


[0,2] 


all 


180 


0.7673 


Pb-Pb 


158 


[-2,2] 


negative 


310 


0.42 


Pb-Pb 


200 


[-0.85,1] 


all 


509 


0.174 


Pb-Pb 


200 




all 


601 


0.1208 


Pb-Pb 


200 


[-1.3,2] 


all 


846 


0.0560 


Pb-Pb 


200 


[-1.7,2] 


all 


1214 


0.0267 


Pb-Pb 


200 


[-2,2] 


all 


1542 


0.01186 



The resulting average multiplicity (N) and entropy index /j,2 are listed in Tables II, III 
and Fig. 4. Also listed in the tables are the colliding nuclei, the incident energy, the particle 
type and the rapidity region used in the analysis. The pt and ip regions in all cases are [0,10] 
and [0, 2ir], respectively. The impact parameter takes a value between and 0.5 fm. 

It can be seen from Fig. 4 that fi2 versus (N) from both Fritiof and Venus Monte Carlo 
simulations fits very well to that expected from the 3-D pure-statistical-fluctuation model, 
independent of the event generator, colliding nuclei, incident energy, particle type and phase 
space region used in the calculation. This means that, in the framework of Fritiof and/or 
Venus event generators, even in the central collision of heavy nuclei at energies up to 200 
A GeV, the statistical fluctuations still dominate the erraticity behaviour. No dynamical 
fluctuation can be observed through erraticity analysis. 

This disappointing fact, however, provides us a possibility to signal the appearance of 
novel physics. The point is that, within the framework of traditional high energy nuclear 
physics the dominance of statistical fluctuations in a given physical process does not depend 
on the concrete conditions, e.g. the collision energy, the mass of colliding nuclei, the cut of 
phase space, etc.. This dominance will disappear and the observed erraticity will deviate 
from that of pure statistical fluctuations only if the events of the studied sample are coming 
from some new kind of physical processes. For illustration, we plot in Fig's. 3 and 4 the 
results from the Gaussian-alpha model proposed in Ref. [12 as stars. It can clearly be seen 
that they do not lie on any of the three curves in these figures. Therefore, we conclude 
that erraticity method has the peculiar property that it is able to filter out all the concrete 
physical conditions used in data analysis and is sensitive to the appearance of novel physics 
in the central collisions of heavy nuclei. 
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